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Abstract 

A new recursive procedure of the calculation of a restricted partition function is suggested. 
An explicit combinatorial formula for the restricted partition function is found based on this 
procedure. 
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1 Introduction 

The problem of partitions of positive integers has long history started from the work of Euler 
who "laid a foundation of the theory of partitions" |2|], introducing the idea of generating 

functions. Many great mathematicians, like Cayley, Sylvester, MacMahon, Ramanujan, and others 

contributed to the development of the theory, using Euler idea. 

Cayley Q found explicit formulas for number pk(n) of partitions of positive integer n into at 

most k parts with small k. He also suggested a method of decomposition of the corresponding 

generating function 

k i oo 

G *(*) = IlT37i = 5>*( n ) 

i=l n=0 

and gave the combinatorial formula for such decompositions (unfortunately, this formula itself 
requires knowledge of all partitions of k). 

Sylvester was the next mathematician who provided a new insight and made a remarkable 
progress in this field. He introduced 0] the so-called Ferrers graphs for presentation of partitions. 
He also found [||, |6| the procedure enabling to determine a restricted partition functions, and 
described symmetry features of such functions. The restricted partition function p(n, d m ) 
pi ii . {d\, d2, ■ ■ ■ , d m }) is a number of partitions of n into positive integers {d\, d%, . . . , d m }, each 
not greater than n. It is very simple to show that the generating function in this case takes the 
form 

m -. oo 

G ( d?n < = II 73^7 = E P(", dm ) tH ■ 

i=l n=0 
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Sylvester showed that the restricted partition function may be presented as a sum of "waves" , each 
wave closely related to prime roots of unit of degree n, where n are prime divisors of elements of the 
set d m . This fact was known to Herschel R who introduced a notion of circulator and Cayley who 
used its elegant version called prime circulator (see || for more information). Namely, Sylvester 
showed that each wave Wj, where i runs over distinct factors in d\, da, ■ d m , is a coefficient of t^ 1 
in the series expansion in ascending powers of t of 

m -. 

^' II > ^ = 2ni P -^+t , u h = 2^-t, (2) 

&ndpi,p2, ■■■iPmax k are integers (unity included) smaller than i and prime to it. It should be noted 
here that the above result is only a recipe for calculation of the partition function and doesn't 
provide an explicit formula. 

Sylvester found || that the shifted partition function 

-, m 

q(n,d m )=p(n--J2di,d m ) 

1 i=i 

has following parity properties: 

q(n,d 2m ) = -q(- n ,d 2m ), q(n,d 2m+1 ) = q(-n,d 2m+1 ), 

and established that these functions have zeros at all integer values of n from to m/2 — 1 for even 
m and at all semiinteger values from 1/2 to m/2 — 1 for odd m. He suggested to use knowledge of 
partition function zeroes for its construction using the method of indeterminate coefficients. 

Recently, a different presentation of the partition function, which may be called polynomial 
expansion, was introduced in ||, where a new recursive procedure for calculation of the restricted 
partition function is found. It permits to reconstruct in an unified way nearly all terms of expansion, 
except one, which also demands usage of the method of indeterminate coefficients. 

In this article I present an approach enabling to overcome this inconsistence, and to determine 
all terms of expansion in the framework of a single recursive procedure. This method also produces 
the explicit formula of the restricted partition function. 



2 Polynomial part of the partition function 

J.J. Sylvester showed that the restricted partition function may be written as a sum of "waves", he 
found the recipe for calculation of each such wave. We consider in this section purely polynomial 
part of the partition function which corresponds to the wave W\. It may be found as a coefficient 
of i -1 in an expansion of the generator 

e st 

GlM) = n£i(i-e-**r (3) 

Sylvester found that W\(s) depends on Bernoulli numbers and sums of powers of elements di. 
We find an explicit form of polynomial part through the Bernoulli polynomials of higher 
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order. We start from the generating function for the Bernoulli polynomials of higher order 
Bt\s\d 1 ,d 2 ,...,d m ) 0: 



St J.71 



where we use a shortcut notation 

d m = {d 1 ,d 2 ,...,d m }. 

One immediately obtains a presentation of W\ generator in the form 



(4) 



1 °° +n 

G 1 (s,t)= - r Y B^\s\ -d m )-, (5) 
v ' ; vr(d m )t m ^ n 'nV 

v ' n=0 



where 

m 

vr(d m ) = n^- 
i=i 

The coefficient of 1 jt in the above expression is given by the term with n = m — 1 
It is useful to consider a shifted restricted partition function defined as follows 

1 m 

V( S ,d m )=W(s-Z(d m ),d m ), £(d m ) = -J>, (7) 

Z i=l 

and its polynomial part Vi(s,d m ) which is cast in 

Vl( S , d™) = {m _i Mdm f inMs - e(d-)| - d™). (8) 

Now we may use formula for Bernoulli polynomials of higher order found by Nor hind [11| 

m 

B^(s\-d m ) = Bir ) (s + Y / d i \d m ' 

=1 



to arrive at 



Then we apply another formula 

n n (m) 
~2< 



i?i m )( g |d m ) = ^Ci 8 j ( S -C(d m ))»^ 



=o 



where £>{ m) (d m ) = 2^ (m) (£(d m )|d m ) and D^ 1 (d OT ) = and obtain 

m— 1 n( m ) 



y '"- d ")= (m -iU-) § ci -' m ""'- L ?- J - (10) 



The quantities £>i m) (d m ) can be calculated using a recursive relation 

n 

D M (cn = c^LA^'V^ 1 ), (11) 
z=o 

where -D; is expressed through the value of Bernoulli polynomial of order I at fixed value of 
argument. 

A = 2^(1/2). 
Instead of one may use more symmetric form 

n m 

^ l) (d m ) = E^Il^.. (12) 

i=l 

where 

n! 

n n^i 

is a multinomial coefficient and summation in (|l2j ) is performed over all such that 

m 

i=l 

The resulting symmetric expression for the polynomial part of the shifted restriction function is 

-i m— 1 I m 

^' dm ) = ^ E CL-i'^E ^11^^(1/2). (13) 



(m - l)!vr(d m 



i=0 i=l 



The inner sum in the above expression can be rewritten using the symbolic notation accepted in 
theory of Bernoulli polynomials (see, for example, |llfl ): 



l m / m \ ' 

£ Of #1^(1/2) = 5>«B(l/2) , (14) 

i=l \i=l / 

where powers rj of l i3(l/2) are converted into indices 

^(1/2)^,(1/2). 

Using this notation one immediately arrives at the compact expression of V±(s, d m ) 

1 / m \ m-1 

tl '^» >-iMd.) l ,t ^ t Ml (15) 



i=l 



Finally, using ([?]), we return to the expression for Wi(s,d r ' 



m—l 



^ dm ) = (m -l M d-) ( a + g^[V2+ fc(l/2)]J • (16) 

It is easily checked that the expression (|D]) verifies the general recursive relation for the shifted 
restriction function (which is valid also for its polynomial part): 

V(s, d m ) - V(s - dm, d m ) = V(s - ^ }d m - 1 ) . (17) 
Introducing the power expansion of V\ in the form 



Fi(s,d m ) = Y,Rf sm ~ j > 



one may write for the coefficients R 



in 



'3 

/ m \ J-l 



R r = (m _^ m A ^nim) • as) 



\i=l 



(m - l)!vr(d'' 

The recursion relation (|Tl| ) is equivalent to the following recursion (1 < j < m): 

R? = — E^cL-i-i+im/m?-! 1 - (19) 

m — j ^ J J 

It should be noted an important distinction between two last expressions - the formula (18) is the 
explicit expression for the polynomial part of the shifted restriction function, while the recursion 
presents an imcomplete procedure (note that j in (|l9| ) cannot be set equal to m, so that R™ is 
not defined in the framework of this procedure) . 

The recursion ( |l9|) is a consequence of ([l7]), which repeated usage leads to more general recur- 
sive relation 

(5m — 1 

V(s + T m ,d m ) = V(s,d m ) + ]T V{s + T m -\ p -d m ,d m ~ l ), (20) 

p=0 

where 

■V = P + 1/2, <5m — T m / d m , 

and 

r m = r(d m ) = LCM(d m ), 

where LCM(d m ) denotes a least common multiple of the set d m . The function V(s,d m ) might 
also be written in a "polynomial" form 

m 

V{s, d m ) = RT(s)s m ~ j , (21) 
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which leads to a recursive formula for r m -periodic function R l J l (s) for 1 < j < m (see 0] 



Rf(s) 



1 



j'-i 



X n d r . 



■ J2 dt l C l m ^_ J+l ■ <f 1 E - =) ■ RJ-i 1 ^ - A P • d m ) ■ (22) 



m — J 7 n 

J 1=0 p=0 

Introducing the shift operator 

S( S) A):{S( S ,A)/( S ) = /( S -A)}, 
we rewrite the above recursion in the form 



1 



E B Ki 

p=0 



m-j — 

Comparison of ( |24| ) with (19) suggests that a replacement 



-)S(s, X p d m j 



(23) 



(24) 



%(l/2) % t = (n/d^- 1 'E ^(l 

7H=0 



(25) 



may be useful in transition from the formulas for the polynomial part of the partition function to 
those of for the function itself. 

Setting in ( ^4|) all BJJ 1 independent of s and using the multiplication theorem for the Bernoulli 
polynomials [10 



m— 1 



E B n (x + -) = m-^B n (mx), 



r=0 



we immediately reproduce (19). The recursion ( p2| ) fails to produce i?™(s), nevertheless partial 
information about it can be extracted. It is useful to separate -R™(s) into two terms 



O0 = 7C(*) + CW, 



where 



TO— l j— i S m — i 

*C0O = E -y- E W 



(26) 
(27) 



2=1 p=0 

is r m _i-periodic function, and the other term, (i m -periodic function, 

r m( s ) = r m( s ~ dm) 

remains indeterminate. In || r™(s) is found by applying the method of indeterminate coefficients 



using knowledge of partition function zeroes. It can be checked using (21) that the partition 
function may be written also as 



m-l 1-1 S m -1 , , ^ , \ 

V(s,d m ) = rZ(s)+T, E B l ( S -±^)lCzl(s + \ p d m ). 

1=1 1 p=0 \ T rn / 



(28) 
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The corresponding expression for the polynomial part of the partition function reads 

^dt 1 _ /l 



(29) 



Setting in the above expression s = we obtain polynomial analog of (p6| 



C = C+E^W:?' (30) 
1=1 1 z 



m 



Substituting into it the general expression (^) we determine r, 

fm-\ \ m ~ 1 

(m-l)!7r(d"») V§ 



T^(£^(V2)J . (31) 



r m 



3 Calculation of r™(s) 

Our goal is to find yet indeterminate r™(s) using ( |3l| ) and recursive relation ( |17| ) for the partition 
function and its polynomial part and its corrolaries ( |l9|) and (p2|) . 

In a simplest case m = 1 it is easily checked that the shifted partition function V(s, {di}) is a 
c?i-periodic function given by 

V(s,{d 1 }) = R\(s)=r{(s) = * dl (s-^)= £ B (^i)^ 1 ( S -A Pl d 1 ), (32) 

z pi=o Tl 

Here Ti = d±, so that pi takes only zero value; we also use parity property of Bernoulli polynomials 
Bq(1 — x) = Bq(x) = 1. The periodic function ^^(s) is defined as a sum of prime roots of unit of 
degree d\\ 

T , , 1 ^ (2iriks\ f 1 , s = (moddi) 

di(s) = ^S exp v^ri = i o . ( m ° drf i) 

The polynomial part of this function 

Vi(s,{d 1 }) = ± 
«i 

this also follows from ( |l5|) for m = 1. 

Consider m = 2, and find Rf (s) using the recursive relation (22) for the set {di, cfe} 



P2=0 



' '^(l-^V 1 ^ ^(i-^^^-A^-A^). (33) 



T2 P2 =o T2 Pl =o r i 
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The symmetry of the problem w.r.t. the permutations of the set elements implies that one can 
apply an interchange d% <-> d 2 in order to produce another valid form of Ri(s). The corresponding 
polynomial part is found as 

p2 _ ij? (l/2) 2 B (l/2) 

It is clear that (^) can be produced from the above as follows - we replace the fraction l/d\ by 
its counterpart *&di( s ~ d\/2) and then use the replacement (|^). The another possible form is 
found by the replacement \/d 2 ^d 2 ( s ~ d 2 /2) and application of (p5|) to it. The free term R 2 (s) 
is constructed in two steps. At first we find TZ 2 (s) assuming d\ the first element of the set. Then 

@3) S ives 



K\{s) = ^f; 1 B 1 (l - ^)Rl(s - X P2 d 2 ) 



P2 =o T2 

T2 £ '^(i-^)' 1 ^ ^(l-^)*^.-^- (34) 

P2=0 T2 pi=0 Tl 

The other part - ^-periodic function r|(s) may be produced in this case by the interchange 
d\ <-> c?2 in the above expression, which corresponds to choice of d 2 as the first element of the set. 

r\{s) = T21 ^ B 1 (l-^)R\(s-\ P2 d 2 ) = 

Pi=0 



pi=o T2 P2 =o Tl 

The corresponding polynomial parts reads 

2 = ^i(l/2) d! 

T2 d x d 2 

and is equal to zero. Nevertheless, we may use the above expression in order to produce (^) - 
we start from the replacement l/d 2 => Vl/<f 2 (s — d 2 /2) and then apply (f25|) to arrive at ^-periodic 
function. 

We conjecture the following formula for (f m -periodic function r™(s) using its polynomial part 
(31) which we rewrite in the expanded form 

i -i m— 1 m 

C = ^ • ^3I)T E II ^ ^(V2). (36) 



i=l 



We start with l/d m =4- ^^(s — d m /2) securing the proper periodicity. Then the replacements 
^5| ) lead to the form 

1 m— 1 m— 1 

CW = 7 TTT E ° r m-i II 'B ri *d m ( S - dm/2). (37) 

(m — 1)! — 
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Using the explicit form of the operator *B r . we arrive at 

TO— 1 m— 1 Ti/di—l 



1 



(m — 1) 



iuE<£-iIK i_1 E ^(i 

r i=l Pi =0 



\ j TO 



(38) 



i=l 



It should be underlined here that the values of "periods" t% depends on the order of the elements in 
the set d m , and in the above formula we have t\ = LCM(d m , d\), T2 = LCM(d m , d±, cfo), ■ ■ ■ , r m _i = 
LCM(d m ),T m = d m , so that p m = 0. The last formula can be rewritten in the symbolic form 

m— 1 



1 



TO— 1 



' i=i 



TO— 1 [Ti/di — l 

E^l E ®(i 

p*=o 



i=l 



(39) 



i=l 



Hence, combination of the last expression with the recursion (25) provides a new procedure for 
calculation of the restricted partition function V(s,d m ). 



4 Explicit formula for restricted partition function 

The same approach can be useful in order to produce a formula for V(s,d m ), starting from its 
polynomial part (13) which can be written as 



TO— 1 



1 lilt, \ lib ± 1 III \ • 

Vi(s,d m ) = (IK 1 ] E CL-iS- 1 - 1 ($> «B(1/2)J . (40) 

It is clear that there exist several equivalent forms of partition function (like two possible forms of 
Ri(s) discussed above), and symmetry considerations can help in selection of the most symmetric 
form. In the discussed case the form is chosen by a selection of the factor 1/di for the replacement 
\/di=> ^fdi( s ~ di/2) and application of (|25|) to the result. 

In general situation we have m choices of 1/di to start with, and because the result doesn't 
depend of such choice, it is natural to seek an expression symmetrized w.r.t. the starting element 
dj. The latter problem is eqivalent to the problem of presentation of the symmetric polynomial 

(TO \ ' l m 

£* = E^IK l 
1=1 / 1=1 

as a sum of m symmetric polynomials, each of them missing only one of terms d{. Introducing 
a function z(r) counting number of zero components of the vector r, one can write the above 
expression as 

TO I -i m 



i=l r Z \ T > n=l 



EE^nc- (4i) 

Using this result we have 

/TO \ / TO 

ik 1 e^^)] =EiEi^n^x(i/2). (42) 



m I m 



\r :J / \,:_] / ! r 
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The polynomial part (flC|) suitable for conversion into V(s, d m ) has the form 



-■ m— 1 m -I I -I m 



(43) 



The explicit expression for the restricted partition function reads 



m— 1 



u( s ,d" 



-| //t _L III 

Z=0 i=l 



71=1 



V^s-di/2). (44) 



Using the actual form of the operators n B rn (25) it is rewritten as 

m— 1 



1 



f)| E Cm-l 5 



i „m— 1— Z 



(m-1). ;=Q 



(45) 



E 

i=i 



r Z y T > n=l pn=0 T ™>* 



n=l 



It should be noted here that the values of "periods" r depends on selected value dj what is reflected 
by an additional subscript i: 



ri,i= di 75,i = LCM(di,d 2 ) 

ri,2 = LCM(di, d 2 ) 72,2 = ^2 

rig = LCM(di, d 3 ) r 2 ,3 = LCM(di, da, d 3 ) 



73,1 =LCM(di, d 2 ,d 3 ) 
r 3 , 2 = LCM(di, d 2 ,d 3 ) 

7"3,3 = d 3 



n, m = LCM(di,d m ) T 2) m = LCM(di, d 2 ,d m ) r 3 , m = LCM(di, d 2 , d 3 , d„ 



r m ,i = LCM(d r ' 
r m , 2 = LCM(d r ' 
r m , 3 = LCM(d r ' 

T~m,m — d m 
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